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Several nonlinear spectroscopy experiments which employ broadband x-ray pulses to probe the coupling 
between localized core and delocalized valence excitation are simulated for the amino acid cysteine at the 
K-edges of oxygen and nitrogen and the K and L-edges of sulfur. We focus on two dimensional (2D) and 
3D signals generated by two- and three-pulse stimulated x-ray Raman spectroscopy (SXRS) with frequency- 
dispersed probe. We show how the four-pulse x-ray signals k\ = —k\ + k 2 + k 3 and fen = fei — fc 2 + fc 3 
can give new 3D insight into the SXRS signals. The coupling between valence- and core-excited states can 
be visualized in three dimensional plots, revealing the origin of the polarizability that controls the simpler 
pump-probe SXRS signals. 



I. INTRODUCTION 

Ultrashort x-ray pulses may be used to probe va- 
lence excitations prepared by a stimulated x-ray Raman 
process. 1-4 In this technique, an inner-shell (core) elec- 
tron is excited to the valence band by an ultrashort laser 
pulse. Then, after the molecule evolves in this highly ex- 
cited state for a time limited by the pulse duration, the 
same pulse stimulates the emission of a photon and the 
core-hole is filled. The final state can be any valence- 
excited state (or the ground state) lying within the pulse 
bandwidth. Because the core to valence transition fre- 
quency is characteristic to the element from which the 
core electron is being excited, this technique can be spa- 
tially selective. For example, in a molecule containing 
a single nitrogen atom, Raman excitation by a pulse 
tuned to the nitrogen core-edge (~ 400 eV) excites only 
those valence-excited states perturbed by the presence of 
a core-hole localized at the nitrogen atom. 

Stimulated Raman is one component of the pump- 
probe signal, which measures the response by exciting 
the system with a short pump pulse, waiting for a pe- 
riod t and then measuring the transmission of the probe 
pulse. The central frequencies of the two pulses con- 
trol which dynamics are probed during the interval, r. 
Shorter pulses and more precise time measurements have 
allowed dynamics at shorter length and time-scales to be 
detected as the technology has improved. 

Time-resolved spectroscopy 5 has been widely employed 
in nuclear magnetic resonance (NMR) to detect tran- 
sitions between nuclear spin levels split by an external 
magnetic field, and is sensitive to dynamics on the /is 
timescale. 6 Spin states have long lifetimes and are easily 
manipulated by external fields. Linear techniques were 
extended to the nonlinear regime, where multiple pulses 
are used to prepare and probe nonstationary states. 6 One 
landmark was the spin echo technique, 7,8 where a simple 
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two pulse sequence is capable of discriminating between 
homogeneous and inhomogeneous line broadening, due 
to fast and slow frequency fluctuations, respectively. 

The photon echo, 9 an extension of the spin echo tech- 
nique into the optical regime, has been applied to molec- 
ular electronic transitions 10,11 and later to vibrational 
systems such as the amide bond stretch in proteins 12,13 
and the hydrogen bond network in water. 14 Pulses with 
central frequencies tuned to optical transitions have been 
used to investigate photosynthetic complexes, 15,16 and 
semiconductors . 17 

The development of x-ray free electron lasers 
(XFELs) 18,19 and higher harmonic generation 20 sources 
may enable pump-probe and photon-echo experiments at 
x-ray frequencies. Zholents and Penn have proposed a 
method to generate attosecond pulse pairs, each mem- 
ber of which is independently tunable through the soft 
x-ray region. 21 Their simulations predict that the tech- 
nique will be able to produce <~ 250 as pulses (bandwidth 
8.5 eV) with center frequencies at the nitrogen and oxy- 
gen K-edges, like those used the simulations presented 
in Section III. Guimaraes and Gel'mukhanov simulated 
infrared (IR) pump/x-ray probe in a model diatomic 
molecule. 22 They proposed that core-hole derealization 
in homonuclear diatomics can be inferred from this sig- 
nal. In a model system with temporally overlapping op- 
tical pump and x-ray probe pulses, it was shown that 
when the probe pulse duration is shorter than one cycle 
of the pump field, the signal is highly dependent on the 
absolute phase of the optical pump. 23 

Optical pump/ x-ray probe experiments have already 
been reported in atoms. 24 Recent work on electron 
and hole mobilities following ultrafast photoionization 
has also stimulated great interest in nonlinear x-ray 
experiments. 25-27 Currently, most studies use the collec- 
tion of molecular ion or photoclectrons to measure inter- 
actions with a probing beam, rather than the change in 
intensity of the probe. 28 Other possibilities which incor- 
porate noisy x-ray pulses have been proposed. 29 

Two-pulse pump-probe experiments, in which the sig- 
nal is defined as the total integrated intensity of the probe 
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with the pump minus that without the pump, were sim- 
ulated in Refs. 1 and 3. These signals were integrated 
over all frequencies, and recorded versus the interpulse 
delay. Interaction with the pump pulse will create core- 
excitations through direct absorption and valence exci- 
tations via a Raman process. By selecting delay times 
which arc longer than the lifetime of the core-excited 
state (< 10 fs), only valence contributions to the signal 
survive. This integrated two-pulse stimulated x-ray Ra- 
man spectroscopy (I2P-SXRS) signal is a direct probe 
of the valence response to a localized core hole which is 
switched on and then off again during the x-ray pulse. A 
three-pulse (3P) extension I3P-SXRS, proposed in Ref. 
2, has two interpulse delays, and thus two dimensions. 
Information on the core-excited intermediate state, the 
higher-lying scattering state, can only be obtained in- 
directly in I2P-SXRS and I3P-SXRS. Increased infor- 
mation is available from these pump-probe signals by 
frequency-dispersing the probe pulse rather than record- 
ing only its integrated intensity, giving an extra dimen- 
sion to the signals and allowing a look inside the probe 
polarizability. 

Frequency-domain x-ray Raman spectroscopy, known 
as resonant x-ray inelastic scattering (RIXS), is a well- 
established technique. Here excitation is achieved us- 
ing a monochromatic beam, and the spontaneous emis- 
sion spectrum is recorded as a function of excitation fre- 
quency. The simplest time-domain x-ray Raman tech- 
nique, I2P-SXRS, uses two pulses, the pump and probe, 
and measures the pump-induced change in the transmit- 
ted intensity of the probe, integrated over all frequen- 
cies. The signal is recorded versus the delay time be- 
tween pump and probe. Due to the vastly shorter life- 
times of core holes compared with valence excitations, 
we can solely focus on valence excitations by examining 
the I2P-SXRS signal for delay times longer than the core 
lifetime (~ 10 fs). I2P-SXRS has additional capabilities 
that RIXS does not. It is possible to use a two-color 
setup where the pump and probe are resonant with dif- 
ferent core transitions, revealing those valence excitations 
perturbed by both core holes. Further, the signal arising 
from the isotropic polarizability can be selected by using 
a magic angle configuration of pulse polarizations. 

Here we add an extra dimension to these techniques 
by measuring the pump-induced change in probe trans- 
mission, as a function of delay time and the probe fre- 
quency. D2P-SXRS shows which core-excited states are 
coupled to which valence-excited states by peering inside 
the probe-pulse effective polarizability. 



II. NONLINEAR X-RAY SIGNALS 

Here we outline the signals considered in this work. In 
four-wave mixing, the molecule interacts with four short 
pulses, and the signal is given by the change in transmis- 
sion of the fourth pulse induced by the other pulses. The 
signals generated in the directions k\ = —k\ + k 2 + fcs 
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FIG. 1. (Left) Pulse sequence for the four-wave mixing setup. 
(Right) The contributing loop diagrams for the fci and fen 
signals. 



and fen = k\ — k 2 -\- k% are represented by the loop di- 
agrams in Fig. 1. In these diagrams, the left and right 
branches of the loop represent the time evolution of the 
ket and bra respectively, and real time flows from bottom 
to top. 30,31 Interactions with the field are represented by 
arrows into (absorption) and out of (emission) the dia- 
gram. The signal is recorded versus the three delay times 
ti, t 2 , and t3, and subsequently Fourier transformed to 
give the 3D signal. Since the system is in a core-excited 
coherence during the ti and t 3 periods, the conjugate fre- 
quencies Oi and f^3 will show resonances in the hundreds 
of eVs. During the t 2 period, known as the waiting time 
or population time, the system is in a valence-excited 
coherence and therefore Sl 2 resonances occur at valence 
frequencies (between 5 and 12 eV for the cysteine model 
used here). Expressions for the k\ and fen signals are 
given in Appendix A. 
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FIG. 2. Loop diagrams for the two-pulse SXRS signal. 

In a two-pulse SXRS pump-probe experiment (Fig. 
2), the system interacts with two pulses with a single 
time delay. The integrated intensity change of the probe 
pulse recorded versus the time delay t 2 gives the one- 
dimensional I2P-SXRS signal. A two-dimensional signal, 
herein called D2P-SXRS, is obtained by further record- 
ing the dispersed spectrum of the probe. The diagrams 
in Fig. 2 are the same diagrams contributing the fci and 
fen signals, but with t\ = 0. Signals are displayed versus 
0,2, the Fourier conjugate of the single time delay. VL 2 
peaks will show valence excitations. Expressions for the 
I2P-SXRS and D2P-SXRS signals are given in Appendix 
B. 

The three-pulse SXRS signals, both dispersed and in- 
tegrated, are represented by the diagrams in Fig. 3. The 
I3P-SXRS signal is defined as the integrated transmitted 
intensity of the probe pulse with two prior pumps, minus 
that without. The two inter-pulse delays, t 2 and £4, con- 
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FIG. 3. Loop diagrams for the three-pulse 2D-SXRS signal. 
SI2 and Q4 are the Fourier conjugates of the two delay times, 
and f2 5 is the monochrometer frequency. Peaks along £^2 and 
CI4 occur at valence frequencies, and those along ils at core- 
excitation frequencies. 



stitute the dimensions of the signal. A third dimension 
can be recovered by frequency dispersing the third pulse. 
Expressions for the I3P-SXRS and D3P-SXRS signals are 
given in Appendix D. 



III. SIMULATIONS 

A. Setup 

We have simulated the spectroscopic techniques shown 
in Figs. 1-3 for the amino acid cysteine. This small 
sulfur-containing molecule provides an important struc- 
tural function by connecting different regions of proteins 
through disulfide bonds. It has been implicated in bi- 
ological charge transfer in the respiratory complex I. 32 
The optimized geometry of cysteine was obtained with 
the Gaussian09 package 33 at the B3LYP 34 < 35 /6-311G** 
level of theory. All restricted excitation windown (REW) 
TDDFT calculations and transition dipole calculations 
for the core excited states were performed with a lo- 
cally modified version of NWChem code 36 ' 37 at the CAM- 
B3LYP 38 /6-311G** level of theory, and with the Tamm- 
Dancoff approximation. 39 For more computational de- 
tails please see Ref. 3. 

Simulations were carried out for four core edges: the 
K-edges of nitrogen, oxygen, and sulfur as well as the 
L2-edge of sulfur. Fig. 4 shows the calculated x-ray ab- 
sorption spectra for these four spectral region. In this 
plot, we have used a Fermi's Golden Rule expression 
convoluted with a Lorcntzian lincshape with a constant 
linewidth (in Ref. 3, an energy-dependent linewidth was 
used to improve agreement with experiment) . Due to the 
presence of two oxygen atoms in cysteine, the density of 
core-excited states near the oxygen K-edge is twice as 
high as those at the K-edges of nitrogen and sulfur, as 
can be seen in the stick spectra of Fig. 4. By examining 
the CI coefficients (not shown here), we see that for any 
given oxygen core-excitation, the core hole is found to 
be entirely localized on one oxygen atom or the other. 
In other words, there is no detectable coupling between 
the two cores at this level of theory. The sulfur L2 edge 
likewise has a higher density of states due to the three 
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FIG. 4. Simulated x-ray absorption (XANES) spectra for 
cysteine at the nitrogen, oxygen, and sulfur K-edges, and the 
sulfur L2-edge. Stick spectra are shown in black, and simu- 
lated spectra with constant Lorentzian linewidths are shown 
in red. The following linewidths were used: Tnis = 0.09 eV, 
Tois = 0.13eV, Tsis = 0.59eV, and r S2p = 0.054eV. We 
plot also the power spectra for the ultrafast laser pulses used 
in the simulations here with FWHM in time and frequency 
of 128 as and 14.2 eV, respectively. The central frequencies 
for the Nls, Ols, Sis, and S2p pulses are 406.4 eV, 532.2 eV, 
2473.5 eV, and 164.6 eV, respectively. 



2p orbitals of sulfur. 

We assume transform-limited Gaussian pulses, with 
full-width half max (FWHM) in intensity of 128 as (14.2 
eV). The power spectra for the pulses are shown as blue 
shaded regions overlaying the XANES spectra in Fig. 4. 
The center frequencies were chosen to cover a significant 
portion of the absorption spectrum, and the large band- 
width can impulsively excite valence-excited states be- 
tween 5 and 12 eV. The calculated UV absorption for 
this model system was given in Ref. 3. 

The coupling strength between valence and core excita- 
tions depends on the dominant hole-particle orbital pairs 
of the excitations. If they share the same dominant hole 
and/or particle orbitals, their coupling is strong and thus 
the corresponding cross peak is visible in the 2D photon 
echo or SXRS spectrum; otherwise the cross peak cannot 
be seen. So by examine the cross peaking pattern in the 
2D photo echo or SXRS spectrum, we are able to gain 
insight about electronic structures of excited states. 



B. Four-Wave Mixing 

The calculated k\ signal for cysteine with all four 
pulses tuned to the nitrogen K-edge, the NNNN sig- 
nal, and polarized parallel to each other (XXXX) is 
shown as a 3D contour plot in Fig. 5. Qi and O3 res- 
onances reveal core-excited states (near 400 eV for ni- 
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FIG. 5. 3D contour plots of |Sfej[, the modulus of the fci photon echo NNNN signal with all four pulses tuned to the nitrogen 
K-edge, and all pulses polarized parallel (XXXX). (Left) The full 3D ki XPE signal. The portion of the signal enclosed in the 
red dashed box (multiplied by 5 here to increase visibility) is responsible for the x-ray Raman resonances. (Right) Enlarged 
spectrum of the region 0.2 > 0, highlighted by the red box in the full spectrum. The walls of the 3D box enclosing the contour 
plot show 2D projections of the 3D data, defined by Eq. 1. 



trogen), while ^2 covers valence-excitations accessed via 
the core-excitations. The most intense peak occurs near 
0,2 = OeV (herein labeled the elastic contribution), cor- 
responding to pathways where the second pulse returns 
the system back to the ground state. One could, in prin- 
ciple, minimize this contribution by redshifting pulse ki 
with respect to pulse k\ such that the second pulse is 
only resonant with emission from core-excited states to 
valence-excited states. However, this is not necessary for 
this system, since the valence excitations are spectrally 
far removed from the elastic contribution (calculated va- 
lence excitation frequencies oj g ' g are between 5.7 eV and 
11.5 eV). 

We focus on the region 5 eV < H2 < 12 eV shown in an 
expanded scale in the right panel of Fig. 5. Also shown 
with the 3D contour plot are 2D projections, displayed 
as the walls to the bounding box. These are defined by 

/oo 
dfii5 fer (n 1 ,o 2 ,n 3 ) 
-00 

S fel (0 1 ,n 2 ,fi 3 )= / dn^fix^fia) (1) 

J5eV 

/oo 
-00 

where the region of integration is chosen to exclude the 
elastic component. 

By varying the phase of the fourth pulse it is possible to 
recover the real and imaginary parts of the time-domain 
signal prior to the application of the numerical Fourier 
transform. Fig. 5 displays the modulus of Eq. A3. In 
Fig. 6 we show the real, imaginary, and modulus of Eq. 1. 
For these and other 2D contour plots presented here, each 
signal is individually normalized to lie between -1 and 1, 



and we use a nonlinear scale 40 to highlight weak features 
in the signal. The real and imaginary signals contain 
additional phase information regarding the molecular re- 
sponse. We look to the k\ signal to explain why certain 
valence states show up in the SXRS signals, and this in- 
formation is more easily read from the modulus signal as 
it does not contain dispersive lineshapes. 

The NNNN k\ signal, shown in Figs. 5 and 6, shows 
that the presence of a nitrogen core hole does not greatly 
effect the valence orbitals. For every core excitation we 
see in f2 l5 there is a single valence excitation to which 
it is strongly coupled. Further, these peaks fall mainly 
along a diagonal line where fii is equal to O2+ ~ 398 
eV. This corresponds roughly to the energy difference 
between the highest occupied molecular orbital (HOMO) 
and the molecular orbital (MO) which corresponds to the 
nitrogen Is core orbital. This is seen most clearly in the 
^3/^2 plot in the bottom row of Fig. 6. Peaks in fl 3 
correspond to the energy difference between core and va- 
lence levels, and we see here that this difference is 398 eV 
independent of ^2- We find roughly a one-to-one corre- 
spondence between valence excitations and nitrogen core 
excitations. We use a picture in which the core electron 
is first promoted to a single unoccupied MO (or a CI 
expansion thereof), following which an electron from the 
HOMO falls down to fill the vacant core hole. A nitrogen 
core hole does not appreciably change the excited parti- 
cle states from those available to valence electrons, and 
the diagonal and horizontal features of the middle and 
lower rows of Fig. 6 show this clearly. 

In Fig. 7 we show the OOOO k\ signal. In this case, 
there is no one-to-one correspondence between valence 
and core excitations. Rather, a given core excitation is 
projected onto many different valence excitations. This 
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FIG. 6. 2D projections of the 3D fci signal of Fig. 5, using an NNNN pulse sequence and XXXX polarization. Top row: 
Sfej (£2i, £^3), Middle row: Sfej(f2i, £^2, ^3), and Bottom row: 5fej(fii, £^2 , ^3). The left, middle, and right columns show 
jSfcJ, KSfej, and 95^, respectively. Each signal is plotted using a nonlinear scale shown in the color bars. 



indicates that a core hole on one of the oxygen atoms 
greatly effects the valence MOs, causing considerable or- 
bital rotation which is not the case for a nitrogen core 
hole. 

The NNNN and OOOO photon echo signals described 
above are single-color experiments. As we can see from 
Fig. 6, the information regarding coupling between the 
core- and valence-excited manifold in the and £^3 di- 
mensions is largely the same. A two-color configuration, 
where two different core orbitals are accessed, allows to 
probe valence-core coupling at different locations. There 
are three possible combinations of four pulses of two dif- 
ferent colors A and B: AABB, ABBA, and ABAB. In Ref. 



41, Schweigert and Mukamel looked at the fei signal from 
different isomers of aminophcnol using the AABB and 
ABAB configuration. Here we focus on the AABB con- 
figuration, as it is the only one that results in a valence- 
excited coherence during the second time interval. Ref. 
41 presents fci signals for constant ti = rather than 
Fourier transforming with respect to ti- 

Fig. 8 shows the 2D projections of the NNOO feu sig- 
nal. The major difference between the fci and fen signals, 
apart from their response to inhomogeneous broadening 
which we do not treat here, is that £^3 resonances occur 
at the core-excitation energies u> eg rather than a,tu> eg >. 
This is due to the fact that the first and second pulse 
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FIG. 7. Same as Fig. 6, using an OOOO pulse sequence and XXXX polarization. 



pair act on the same side of the density matrix, which 
makes the fen signal more straightforward. The bottom 
row of Fig. 8 shows jS^^!, fi 2 , ^3)], an d makes it clear 
that the coupling between the valence-excited and oxy- 
gen core-excited manifolds is much different than that 
for a nitrogen core excitation. The first few oxygen core- 
excited states, with energies between 532 eV and 535 eV, 
are each coupled to many different valence excitations. 
The |Sfe,(fii, f2 2 , f^3)| plot in the middle row of Fig. 8 
is similar to the same plot from Fig. 6. However there 
are distinct differences. The valence-excited state at 5.74 
eV is absent from the NNNN signal, but is strong in the 
NNOO signal. Furthermore, in the NNOO fcn signal, we 
see that this state is coupled to multiple nitrogen core- 



excited states, in contrast to NNNN signal where there 
was a one-to-one correspondence between the two mani- 
folds. 

The projections defined in Eq. 1, and displayed in Figs. 
6 and 8 provide only one example of how the full 3D 
spectrum can be parsed into 2D spectra, which are more 
easily interpreted. Another method, one which will give 
direct insight into the origin of the two-color I2P-SXRS 
resonances, is to display 2D fii/f^ plots for constant il 2 . 
In Fig. 9 we show slices of the OOSS fcn signal for con- 
stant 1^2 corresponding to three different peaks in the 
OS I2P-SXRS spectrum. These give a correlation plot 
between core-excited states located at different atomic 
centers, corresponding to a specific valence-excited state. 
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FIG. 8. 2D projections of the 3D fen signal using an NNOO pulse sequence with XXXX polarization. 



We see that the valence-excitations at 6.6 eV and 8.9 eV 
have are coupled to the same sulfur Is core-excited states, 
but different oxygen core-excited states. The higher- 
energy 11.4 eV valence excitation, on the other hand, has 
a much different excitation pattern, being mainly coupled 
to higher-lying sulfur excited-states. 

Finally we note that the four-wave mixing signals can 
be used to disentangle the effects of pulse polarization 
in x-ray Raman signals. Polarization is more important 
in x-ray Raman than in traditional optical Raman. In 
a vibrational resonant Raman experiment, the system 
is transiently promoted to an electronically-excited state 
before de-excitation to a vibrationally excited state in 
the ground electronic state. Under the Condon approxi- 



mation, the transition dipoles for the upward and down- 
ward transitions will always be parallel to each other. 
In x-ray Raman spectroscopy, there is no constraint on 
the angle between the upward V eg and downward V eg ' 
transition dipoles. In general, for a third-order nonlinear 
experiment such as the fci and fen, there is an orienta- 
tional factor , 2 'v'fV V , multiplying each term in 
Eqs. A3 and A4, defined in Ref. 42. By setting the 
angles 

(ei, e 2 , e 3 , e 4 ) = (0°, 0°, 6 MA , 6 MA ) (2) 

where V is vertical and 9ma — 54.7° is the magic an- 
gle, the orientational factor becomes equal to (V e > g > ■ 
Vge'){V eg i ■ V ge ). This is one of the three intrinsic 
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isotropic signal components. This polarization configura- 
tion corresponds to the 2P-SXRS signal when the pump 
and probe are polarized at the magic angle to each other. 
It is possible to isolate all three of the isotropic compo- 
nents. When 



(ei,e 2 ,e 3 ,e 4 ) = (0°, 9 MA , 0°, M a), 



(3) 



the orientational factor becomes (V e ' g > ■ V e g')(Vg e ' ■ V ge ) 
And finally, when 



(ei, e 2 , e 3 , e 4 ) = (0°, 6> M a, #ma, 0°) 



(4) 



the orientational factor becomes (V e 'g' • V g e){V eg < ■ V ge >). 
The signal resulting from an arbitrary pulse polarization 
configuration can be written as a linear combination of 
these three isotropically averaged components. Fig. 10 
depicts the fen signal, using an OOSS pulse sequence, for 
the three field polarization configurations described here. 

These polarized signals show distinct differences re- 
lated to the angles between the various transition 
dipoles involved. For example, the largest peak 
in the (0°, #ma, 0°, 9ma) signal, shown in the mid- 
dle column of Fig. 10, occurs at (f2i,fi 2 ,f2 3 ) = 
(532 eV, 6.6 eV, 2474 eV) and is mostly absent from the 
(0°, 0°, 0mAj $ma) signal, shown in the left column. This 



indicates that the transition dipoles connecting the 
ground state to the given core-excited states are nearly 
perpendicular to the transition dipoles between these 
same core states and the valence-excited state at 6.6 eV. 
The polarization dependence of broadband x-ray signals 
can be used as an experimental check on the inter-dipole 
angles predicted by electronic structure calculations. 
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FIG. 10. fen signals simulated for an OOSS pulse sequence for three different polarization configurations. Each column 
represents a different polarization schemes, with the polarization vectors of pulses 1 through 4 indicated at the top of each 
column. The top row shows |Sfe n (fir, O2, | , the middle row |5fe n (f2i, fi2, fis)|, an d the bottom row [/fjfe n (fir, Q2, fis)|- 



C. Dispersed versus Integrated Two-Pulse SXRS 

The frequency-dispersed two-pulse SXRS signal reveals 
information on the coupling between valence excitations, 
excited impulsively by the pump pulse, and the core- 
excitations resonant with the probe pulse. The integrated 
SXRS signals for cysteine were presented in Ref. 3. 

In Fig. 11 we show the D2P-SXRS signals, using all- 
parallel polarization, using a probe pulse tuned to the 
nitrogen K-edge and pump pulses tuned to the four core 
edges considered here. We can see from Eq. B4 that any 
given resonance between a valence-excited state g' and 
a core-excited state e, will generate two peaks located 
at (ri 2 , ^3) = {w g 'g,u eg ) and (cj g > g ,oj eg /). Because our 
model does not possess valence excitations with energies 



below 5.74 eV, it is easy to distinguish between these 
types of peaks. In Figs. 12 through 14 we show the re- 
maining D2P-SXRS signals for all other probe detunings. 
Along with the 2D contour plots, we also show the ID 
projections along each axis, defined analogously to Eq. 
1. The topmost projection is the integrated 2P-SXRS 
signal. 

In Fig. 15 we compare the fei and k\\ signals with the 
dispersed Raman signals, both taken with a two-color 
setup. If the J7i axis is integrated out, formally equivalent 
to setting t\ = 0, then the pump-probe signal gives iden- 
tical information to the experimentally more challenging 
fei and k\\ signals. However, in systems which feature 
both low- and high-energy valence excitations, where the 
Lu eg and Lu eg i peaks would overlap, the k\ and k\\ sig- 




6 8 10 12 6 8 10 12 

n 2 (eV) 0.2 (eV) 



FIG. 11. Frequency-dispersed SXRS spectra (Eq. B4)for cys- 
teine, for various pump pulse detunings with the probe tuned 
to the Nls-edge transition frequency. The vertical axis is the 
dispersed frequency and the horizontal axis is the Fourier 
conjugate of the interpulse delay ■ The 2D signal is plot- 
ted using an arcsinh nonlinear scale, and the projections are 
shown, in a linear scale, in the top and right marginals. The 
trace along the top corresponds to the integrated SXRS sig- 
nal. 

nals would be able to more clearly identify resonances 
between core- and valence-excited states. Furthermore, 
four-wave mixing gives direct cross-correlations between 
core-excitations on different atoms while this information 
must be inferred from the pump-probe measurements. 
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6 8 10 12 6 8 10 12 



n 2 (eV) fl 2 (eV) 

FIG. 12. Same as Fig. 11, but with the probe tuned to the 
oxygen K-edge. 
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FIG. 13. Same as Fig. 11, but with the probe tuned to the 
sulfur K-edge. 
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6 8 10 12 6 8 10 12 

n 2 (eV) 0.2 (eV) 



FIG. 14. Same as Fig. 11, but with the probe tuned to the 
sulfur L-edge. 



OONN 4- Wave Mixing ON FD2P-SXRS 




fl 2 (eV) fl 2 (eV) 

FIG. 15. (Left column) Four-wave mixing signals using 
the OONN pulse sequence (XXXX polarized). On top is 
|Sfc n (fii, £12, £^3) | and on the bottom is \Ski (fii, ^2, Os)|. 
(Right column) The modulus D2P-SXRS signal with an ON 
pulse sequence, XX polarized. 



D. Dispersed versus Integrated Three-Pulse SXRS 

In Fig. 16 we show four examples of the I3P-SXRS sig- 
nals for cysteine, with the first pulse resonant with the 
sulfur K-edge and the third pulse at the oxygen K-edge. 
The signals differ in the tuning of the second pulse. Di- 
agonal contributions to the signal, where = ^4, are 
largely insensitive to the identity of the middle pulse, 
while the trace along = are insensitive to the probe 
pulse. That is, for an integrated 3-pulse SXRS signal 
with an ABC pulse configuration (where A,B, and C rep- 
resent different cores to be excited), the diagonal trace 
will be similar to the I2P-SXRS signal with an AC pulse 
configuration and the SI4 = trace equal to the I2P- 
SXRS signal with an AB pulse sequence. In Fig. 17 we 
show the SOO D3P-SXRS signal. 




Ch (eV) n 2 (eV) 

FIG. 16. I3P-SXRS signals from cysteine with the first pulse 
resonant with the Sis transition and the third at the Ols 
transition, with XXX polarization. 

As can be seen in Fig. 3, the three-pulse Raman signals 
involve core-excited states e, e', and e". e and e' con- 
tribute to ctg/g and a gg " (Eq. Dl), respectively, whereas 
e" contributes to a g " g ' (Eq. Dl). In our current approx- 
imation, core excited states are given by |e) = c\c n \g) 
where and c are the Fermi creation and annihilation 
operators, and a and n are virtual and core orbitals, re- 
spectively. While this is a good representation for e and 
e', additional contributions |e") = (? h S a CiC n \g) , where i 
and b are valence and virtual orbitals, respectively, should 
be taken into account for e" . These double excitations 
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fi 3 (eV) 



H 2 (eV) 



FIG. 17. D3P-SXRS signals from cysteine using an SOO pulse 
configuration, with XXX polarization. Q2 and H4, are valence 
excitations, and Q5 represents core excitation. 

are neglected in the present simulation, and should result 
in additional peaks. 

IV. CONCLUSIONS 

We have shown how the fci and fen x-ray four-wave 
mixing techniques can reveal the coupling between core- 
and valence-excited states. Furthermore, it is possible to 
isolate the three independent rotationally averaged con- 
tributions to the signal in these experiments. 

The frequency-dispersed two-pulse SXRS signal, a 
much simpler experiment, gives the same information as 
the photon-echo provided the first time delay is set to 



zero. The D2P-SXRS allows to look inside the probe po- 
larizability, while the four-wave mixing reveals both the 
pump and probe polarizability. 

The time-domain Raman signals have have a fur- 
ther advantage over their frequency-domain counterpart, 
RIXS, due to their directionality. In a RIXS experiment 
the spontaneous signal is not phase-matched and will be 
emitted evenly over a An solid angle, only a fraction of 
which will be collected. In an SXRS experiment the sig- 
nal is emitted exclusively in the probe direction, so the 
entire signal is recovered. This advantage applies to the 
fei and fen signals as well, where the signal is emitted in 
a background-free direction. 

We have also proposed extending the three-pulse SXRS 
experiment by adding frequency-dispersed detection. An 
example of this three-dimensional fifth-order signal was 
given. However, the REW-TDDFT method used to cal- 
culate core-excited states, while sufficient for the third- 
order echo and pump-probe, cannot represent states with 
both core- and valance-holes. The development of quan- 
tum chemical methodologies more suited to this experi- 
ment is called for. 

The XFEL light sources are very bright (e.g. 10 12 pho- 
tons in 85 attoseconds at frequency 8 keV); 43 the first 
experiments performed at the LCLS saturated the core 
x-ray excitation in a molecular beam of nitrogen gas. 44 
This shows that the stimulated Raman process at these 
intensities can compete with the Auger process which 
may be detected. The strength of the nonlinear signal 
has been estimated in Ref. 45 by treating the field semi- 
classically. 
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Appendix A: Four-pulse Phase-Matched k\ and fen Signals 

In this four-photon stimulated process (Fig. 1), the molecule interacts with four short pulses, described by 

E(r,t) = ^2e j £ 3 {t- Tj )e l ^- r + c.c. (Al) 

j 

where ei,£j, kj and Tj are the pulse polarization vector, temporal envelope, wavevector, and arrival time, respectively. 
The signal is given by the change in transmission of the fe 4 pulse induced by the other pulses, and is recorded versus 
the three interpulse delays, defined by t\ = t-i — t\, etc. (see Fig.l). 
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By placing the detector in a particular wave-vector matching direction, it is possible to selectively probe a subset 
of the pathways which contribute to the overall molecular response function. In this paper, we look at the fci and fen 
signals, where we set k 4 equal to k\ = —k\ + k 2 + k 3 and fcn = fci — k 2 + k 3l respectively. The fen signal can be 
obtained by using the same experimental geometry as in the fci signal, if the order of pulses 1 and 2 are interchanged. 

We calculate the fci and fcn signals using time-dependent perturbation theory with the interaction Hamiltonian 

H int (t) = -V ■ E{t) (A2) 

where V is the transition dipole operator. For all of the transitions considered here, the x-ray wavelength is more 
than an order of magnitude larger than the core orbital, so the dipole approximation is justified. 

Contributions to these signals come in three varieties, the ground-state bleach (GSB), excited-state emission (ESE), 
and excited-state absorption (ESA). In the ESE and ESA contributions, the first two pulses act on different sides of 
the density matrix creating a core-excited population which is then probed via stimulated emission or absorption of 
the second pulse pair, respectively. In the GSB term, the first two pulses act on the same side of the density matrix 
in an up-down fashion, creating a valence-excited wavepacket which is then probed by the third and fourth pulses. 
Since the lifetimes of core-excited states are very short (less than 10 fs for the cores considered here 46 ) than those for 
valence excitations, the ESE and ESA terms can be eliminated by restricting t 2 to be longer than the core lifetime. 
We assume well-separated pulses. 

The four-wave mixing signals are given by 

<5fci(— 1; — fi 2 , 3 ) = 

\p £i{u e <g)£2{u e >g>)£Z(u eg >)£ 3 (uj eg )V ge V eg >Vg< e <V e > g (A3) 
~^ , (-fii +u e >g +ir e /)(-Q 2 + ijJ g > g + iT gl )(n 3 -uj egl + iT e )' 

where £i(co) is the spectral envelope of the i th pulse, oj rs is the transition frequency between states r and s, and T e is 
the inverse lifetime. The fcn signal is likewise given by 

Ski{— fii, — ^2, 03) = 

£l(u]e'g)£2(Ue>g')£t(Ueg')£3(Ueg)VgeVeg'Vg>e>Ve'g (A4) 



E 



■f t (-ill + U} e >g + ir e /)(-Q 2 + LUglg + lT g') (Cl 3 ~ CO £ g' + lT £ ) 

Four-wave mixing is a natural extension that sheds light on the SXRS process. The two-pulse and three-pulse 
versions of SXRS arc pump-probe experiments wherein the system interacts twice with each pulse. In the photon 
echo experiment, the pump pulse is split into two, pulses fci and fc 2 , each of which interact with the system just once. 
Likewise the probe pulse is split into pulses fc 3 and fc 4 . 

In the simulations presented in Sec. Ill, we perform ensemble averaging over the random distribution of molecular 
orientations using the tensor formalism. 42 ' 47 



Appendix B: Two-pulse SXRS Stimulated X Ray Raman Signals 

The integrated SXRS signal (Fig. 2) is one dimensional. The signal is given by 1 ' 3 

/oo 
dt 3 £^t 3 )P^\t 2 ,t 3 ) 
-OO 

= -E 



a 'igg' a Ug'g | ( a 2;gg' a i;g'g)* \ l51 ' 

^2 - LOg'g + iFg' ^2 + LOg'g + iT 



where 

otj = a/ + ioij 



e,g',g" M 

is the effective isotropic polarizability averaged over the spectral envelope of the j th ultrashort pulse, £j. The two 
terms in Eq. Bl have peaks in the positive and negative fi 2 regions, respectively, which carry the same information. 
We thus only plot the positive f2 2 region, retaining only the first term in Eq. Bl. 
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Note that the signal depends only on the anti-Hcrmitian part of a 2 - This is because it represents the change in 
the intensity of the probe pulse integrated over all frequencies due to the presence of the pump. If the probe is off 
resonance, then a 2 ' = and its net absorption by the system will vanish regardless of the pump. The transmission 
of an off-resonant probe is a parametric process, energy can be exchanged between different modes within the probe 
bandwidth, but will not change its integrated value. 

A two-dimensional signal can be obtained with two pulses by sending the probe through a spectrometer, and 
recording the dispersed spectrum as a function of the interpulsc delay, which is then Fourier transformed numerically. 
The D2P-SXRS signal is given by 

SD2P-sxRs{t2,n 3 ) = z£Z(n 3 )P {3) (t 2 ,n 3 ) (B3) 

which, after a Fourier transform with respect to the delay time, evaluates to 

o fO O \ \ " iV geVeg'(ai;g'g) ( £-2 (gg + ^ gg' )gj (gg) ^fflg+^Vgj^jg) \ 

^ Q 2 ~ UJg'g + iTgf \ fi 3 - UJ eg + lT e il 3 - LU eg < - lY e J 

It can easily be verified that 



Sl2P-SXRs{^2) = J dil 3 SD2P-SXRs{Q2,Q 3 )- (B5) 

Eq. B4 shows that while the frequency-dispersed signal becomes weaker when the probe is off resonance, it does not 
vanish. However, the integrated signal does vanish, as is shown in Appendix C. 

Appendix C: The Effective Polarizability 

The effective polarizability is defined as 

. — . v , 1/ r°° _ „ 

duj jy ' Jy y y/ (CI) 



_ \- Vg'eV eg f°° A {U) £ * {LU + ^ 



This operator is non-Hermitian, and we may write the matrix elements of its conjugate transpose as 

a )\g'g = ( a r,99'Y 

Vg'eVeg f°° , £j (u) (u + Ugg>) 



^ 2tt uj-uj eg -iT e (C2) 



2tt J_ 00 uj - uj eg > - iT e 

It is useful to explicitly split aj into Hermitian and anti-Hermitian parts as 

a'j =(aj + otj)/1 
ia'j —{oij — ctj)/2 

whose matrix elements are 



(C3) 



^r,g g 



v r due* m e j ( W + w g ,g) ( — ^— + — 1 — —) 

^ Air J_ 00 y y ' \uj -uj eg , +iT e uj - uj eg , - iT e J 



(C4) 



and 



ia j;g'g 



= E^ £ r du£* (.)£,(. + tOg,g) 



V > V f°° / 1 1 

, UJ - UJ eg < +iT e UJ - UJ egl - lT ( 



(C5) 



Ti + {uj-uj eg ,y 
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In the off-resonant case, where £j(uj eg i) = 0, the Hermitian part of oij does not go to zero while the anti-Hermitian 
part does. This can be seen explicitly by taking considering Eqs. C4 and C5 in the limit where T e — > 0. 
We can also easily show that the D2P-SXRS signal vanishes when the probe is off resonance, 



iVgeV e g'{a 1;glg ) 



. --2 — Uqig + iV„i 
e.g ^ ^ u 



^ £ 2 (n 3 + L, gg ,)£S(n 3 ) /■ dn ^(n 3 + u g , g )£ 2 (n 3 ) 



SsxRs(n 2 ) = Y l ^ seVeAa ^ s '" ) [ <m 3 £ 2 (n 3 + uj gg ,)£Z(n 3 ) 

'—'112 — W q i„ +IL „' J 



- J g'g "I" '<"■ g' J V^3 ~ ^eg + i^e ^3 - ^eg ~ ^ c 

^ Vg e V e g,(a l g , ) f dQ £ ^ ) ^ 

Q 2 - W fl ' fl + lT g , J (fi 3 - W eff ) 2 + Tl 



(C6) 



^3 - w e 9 + «r e 7 Q 3 - LU eg < - lY e 

First we make the substitution in the second integral f} 3 = f2 3 + oj g > g , giving 

^,0.2- Ug'g + lYg, \J ^ 3 - W e g + lY e J fig - LU eg - lY e J 

we then combine the two integrals, renaming Q 3 to Q 3 , giving 



(C8) 



The fraction in the integrand will be non-negligible only near resonance, but if the probe-pulse spectral envelope 
vanishes near resonance, the integral will vanish as well. 



Appendix D: Three-pulse SXRS Stimulated X-Ray Raman Signals 



The I3P-SXRS signal (Fig. 3), written in a consistent notation to that above, is 

<? (O n\-V iai '9'a ( a l,gg" a 3;g"g' a 3; gg " a ^g"g' \ mu 

gT^„ "2 - + U g' \ "4 ~ Wg'g" + U 9 ' "4 ~ ^ g "g + U g" y 

Eq. Dl contains only positive resonances. Note the symmetry, Si 3 p-sxrs(~ &2, — ^4) = (Si3p-sxrs{&2, ^4))*- 
To obtain the expression for the frequency-dispersed three-pulse SXRS (D3P-SXRS) signal, we simply need to 
"unpack" the probe polarizability in Eq. Dl, just as Eq. B4 was obtained from Eq. Bl. 



Sd3P-SXRs(^2, ^5) — ^ 



e,g',g" 



l ifl ' fl 4 3 ,g" /£ 3 (Q 5 +W s < V )£ 3 *(ft 5 ) ^(Og+uy^sC^) 



— "I" *Fg' — ^g'g" T *Fg' \ O5 — <-^eg" "F *F e $25 — W e <j' — *F e 

iQ!l ;s ' g a 2;g",g' {£3(^5 + W S9")^3 (^5) £3 (^5 + U g n g )£ 3 (fit 



(D2) 



^2 - ^q'q + TV fi 4 - + iTo" V ^5 — ^eo + i^e ^5 ~ W™" - lY e 
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